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In this letter, we show how the combination of both a gauged U(l)x family symmetry and 
an extended vertical gauge symmetry in a single model, allows for the presence of additional 
Clebsch texture zeroes in the fermion mass matrices. This, leads to new structures for the 
textures, with increased predictivity, as compared to schemes with enhanced family symmetries 
only. We illustrate these ideas in the context of the Pati-Salam gauge group SU (4) <g> SU (2) l <8> 
SU(2) R supplemented by a U(l)x gauged family symmetry. In the case of symmetric mass 
matrices, two of the solutions of Ramond, Roberts and Ross that may not be obtained by family 
symmetries only, are accurately reproduced. For non-symmetric textures, new structures arise 
in models of this type. To distinguish between the solutions in this latter case, we performed a 
numerical fit to the charged fermion mass and mixing data. The best solution we found allows 
a fit with a total x 2 of 0.39, for three degrees of freedom. 



In the recent years, there has been a lot of effort in trying to understand the pattern of 
quark and lepton masses and mixing angles. In order to explain the observed hierarchies in 
the most predictive way, it has been proposed that zero textures in the Yukawa matrices exist, 
like in the Fritzsch ansatz |1] and the Georgi-Jarlskog (GJ) texture 0. Ramond, Roberts and 
Ross (RRR) || have made a survey of possible symmetric textures which are both consistent 
with data and involve the maximum number of texture zeroes. Such textures can arise due to 
a family symmetry group G and some new heavy matter of mass M which transforms under 
G, resulting in effective non-renormalisable operators that generate fermion masses [|J. Such 
family symmetries arise in most of the superstring models. A realisation of this picture has been 
provided by Ibanez and Ross (IR) ||, based on the MSSM extended by a gauged family U(l)x 
symmetry with 9 and 9 singlet fields with opposite X charges, plus new heavy Higgs fields 
in vector representations. Anomaly cancellation occurs via a Green-Schwarz-Witten (GSW) 
mechanism, and the U(l)x symmetry is broken not far below the string scale, generating 
Yukawa matrices of the form 
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These matrices resemble RRR solution 2 || and for e = e 2 , e = 0.22, they reproduce the 
known fermion mass hierarchy. However, the 23 = 32 element in the down mass matrix tends 
to be large, predicting a value for \V c b\ that has to be lowered by the insertion of coefficients 
or by a cancellation due to phases. On the other hand, one can imagine a different fit, that 
would naturally lead to a smaller 23 element. This can be achieved by introducing a small 
parameter 5 which originates from some flavour independent physics and appears as a factor 
in all non-renormalisable elements, so that e.g. the down quark mass matrix is modified to 

/ 5e 8 5e 3 5e 4 \ 
A D = 5e 3 5e 2 5e . (2) 
V 5e 4 5e 1 / 

For e ~ 0.22 and 5 ~ 0.2, in the down quark mass matrix, one finds that all entries have 
naturally the correct magnitude to reproduce the desired phenomenology. 

However, in order to naturally obtain such a contribution 5 in addition to the expansion 
parameter e in the mass matrix entries, one has to go beyond the MSSM. The basic new idea 
of our approach is to combine the ideas of gauged U(l)x family symmetry with unification in 
the form of an extended vertical gauge symmetry. As a concrete realisation of this idea we 
consider a specific supersymmetric unified theory that can be derived from a superstring 
model 0, based on the Pati-Salam gauge group SU(4) ® SU(2) L ® SU(2) K [g]. We emphasise 
that the results can be generalised to different groups; the important point is that the quarks 
and leptons are unified into a common representation, leading to new Clebsch relations. In 
this particular model, the left-handed quarks and leptons are accommodated in the following 
representations (see || for details) 

F* aa =(4,2,l) ,F^ = (4,1,2) (3) 



1 



where a = 1, . . . , 4 is an SU(4) index, a,x = 1,2 are SXJ(2) LjR indices, and i — 1, 2, 3 is a family 
index. The Higgs fields are contained in h x a = (1, 2, 2) while the two heavy Higgs representations 
are 

H ah = (4,1,2), H ax = (4,1,2) (4) 

The Higgs fields are assumed to develop VEVs < H >=< u H >~ M GUT , < H >=< u H >~ 
M GUT , leading to the symmetry breaking at M GUT , SU(4) ® SU(2) £ ® SU(2) R — ► SU(3) C <g> 
SU(2) L ® U(l) y in the usual notation. Under this symmetry breaking, the bidoublet Higgs 
field h splits into two Higgs doublets hi, h 2 whose neutral components subsequently develop 
weak scale VEVs, < h\ >= Vi, < h® >= t>2 with tan/3 = v%/vi. The relevant part of the 
superpotential involving matter superfields is 

W = X^FtFjh + X^FjHOk + \ihh -\ (5) 

where 0j are the superfields associated with the singlets. From Eq.|^ we find that the Yukawa 
couplings satisfy the boundary conditions 

Xi(Mgut) = \%(M GUT ) = \%(M GUT ) = \%(M GUT ) = \%{M GUT ), (6) 

When Eq.(||) is applied to the third family , successful predictions for top, bottom and tau 
masses are obtained. However these relations do not hold for the lighter families. The standard 
way forward is to suppose that there is some family symmetry which prohibits the renormalis- 
able terms for the lighter families, but permits them for the third family. The lighter fermion 
masses are then accounted for by non-renormalisable operators whose order is controlled by the 
family symmetry. 

The non-renormalisable operators which were studied in ref. [H], are formed from different 
group theoretical contractions of the indices in 

= F aa F j3x h y a H lz H p ' w . (7) 

and are of the form: 

On - (FiF^h (j^L 

where M > M G ut (in string models it can be associated with the string scale). When H, H 
develop their VEVs, such operators will become effective Yukawa couplings of the form FFh 
with a small coefficient of order M GUT /M 2 . These operators were previously written down 
independently of any particular family symmetry. Now we extend this scheme by introducing 
a U(l)x family symmetry which is broken at a scale Mx > M G ut by the VEVs of the singlet 
fields 6 and 6. To get an idea of the implications, for simplicity we assume that (i) the singlet 
fields do not couple directly to the matter fields F and (ii) the the combination HH has a 
zero quantum number under the symmetry. The new operators are then modified from those 
in Eq||: 

(H ' Ff\ ( f) n f) m \ 
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where again M' represents an energy scale which can be associated with the string scale or 
the U(l)x breaking scale. The single power of (HH) is present in every entry of the matrix 
and plays the role of the factor 5 in Eq.^], while higher powers of (HH) are suppressed. In 
contrast to the MSSM case, here the factor {HH) carries important group theoretical Clebsch 
information. In fact Eq.|9] amounts to assuming a sort of factorisation of the operators with 
the family hierarchies being completely controlled by the 9, 9 fields as in IR. Here, m and n 
are dependent on i, j, while the splittings between different charge sectors of the same family 
are controlled by the Clebsch factors in (HH). The Clebsch factors have of course a family 
dependence, i.e. they depend on i, j. The relevant n = 1 operators, appear in Table 1. 

We should comment on the origin of the operators in Eq|| There are two possible sources 
of these operators which we can envisage: (1) the operators may emerge directly from the 
superstring construction; (2) they may be generated from the effective field theory below the 
string scale. In case (1) the operators will be suppressed by powers of the string scale, while 
in (2) the operators will be suppressed by powers of the mass of some new heavy states in 
vector-like representations which mix with the quark and lepton representations. The first 
possibility will be addressed in fLOf . Here, we would only like to remark that, in this case, the 
larger contributions naturally arise by contracting fields belonging to the same sector of the 
superstring (Neveu-Schwarz, or Ramond sector), and that the operators that arise in this way, 
lead to correct phenomenological predictions. As for the second possibility, this would require 
heavy vector-like representations of matter fields ip^ x ' + ip(~ x ^ where the X charges are shown 
in parentheses. Flavour mixing occurs via insertions of 9 and 9 fields leading to the so called 



spaghetti diagrams recently discussed in the context of the MSSM In the case of the 



Pati-Salam model, additional heavy Higgs fields in adjoint representations may generate the 
operators in Eq.|| An example of a spaghetti diagram which can yield a 23 entry in the Yukawa 
matrix is shown in Fig.l. Note that the HH pair will be in the Pati-Salam representation of the 
heavy £ field which does not change flavour since it has zero X charge (in parentheses). Flavour 
changing occurs only in the heavy sector (where these fields have the same Pati-Salam 
representations as the chiral matter fields Fi, Fj) via insertions of the 9 fields. 

Note that Table 1 includes cases of zero Clebsch coefficients, where the contribution to the 
up-type matrix, for example, is precisely zero. Similarly there are zero Clebsch coefficients for 
the down- type quarks (and charged leptons). The existence of such zero Clebsch coefficients 
enables us to obtain textures with different places for zeroes in the up and down mass ma- 
trices, something that is not easy in the case that only a flavour symmetry is included in the 
model. Apart from the zero Clebsch coefficients another advantage of combining quark-lepton 
unification with a family symmetry is that one can also account for the mass splitting within 
a particular family, since now the same expansion parameter appears in both the up and down 
quark mass matrices. 

To see how all this occurs, let us take a specific example. Consider the symmetric n = 1 
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Figure 1: A Pati-Salam spaghetti diagram. 
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Table 1: The n = 1 operators with Clebsch coefficients as shown. 
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where O33 is the renormalisable operator and s.d. stands for a sub-dominant operator with a 
suppression factor compared to the other dominant operator in the same entry. Substituting 
the Clebsch coefficients from Table 1 we arrive at the following \ U > D > E Yukawa matrices, at the 
GUT scale, 
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where A^ and A^ arise from the dominant O^, operator and A^ comes from a sub-dominant 
operator that is relevant because of the texture zero Clebsch in the up sector of O^,. The 
zeroes in the matrices correspond to those of the RRR solution 5, but of course in our case 
they arise from the Clebsch zeroes rather than from a family symmetry. The numerical values 
corresponding to RRR solution 5 are, 
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In our model the hierarchy A 22 << ^22 ^ s explained by a Clebsch zero and a suppression factor 
of the sub-dominant operator. Using Eq.[l2| we can read off the values of the couplings which 
roughly correspond to a unified matrix of effective (and dominant) couplings: 




3 x KT 3 1 x KT 3 
3 x KT 3 2 x KT 2 2 x KT 2 
1 x KT 3 2 x KT 2 1 



(13) 



where we have extracted the Clebsch factors. Thus, the Clebsch factors reproduce the values 
in Eq.[T2| required by phenomenology for the dominant effective couplings displayed in Eq.[T3|. 

The case we are examining is different from the IR analysis in two aspects: (a) the fermion 
mass matrices of different charge sectors have the same origin, and thus the same expansion 
parameter and (b) all variations between these sectors arise from Clebsch factors. The structure 
of the mass matrices is again determined by a family symmetry U(l)x, but now the charge 
assignment of the various states are as in Table 0. 
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Table 2: U(l)x charges for symmetric textures. 
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The need to preserve SU(2) L invariance requires left-handed up and down quarks (leptons) 
to have the same charge. This, plus the additional requirement of symmetric matrices, indi- 
cates that all quarks (leptons) of the same z-th generation transform with the same charge aij. 
Finally, lepton-quark unification under S77(4) (g> SU{2) L ® SU{2) R indicates that quarks and 
leptons of the same family have the same charge. The full anomaly free Abelian group involves 
an additional family independent component, U(1)fi, and with this freedom U(l)x is made 
traceless without any loss of generality, giving a.\ = — (a 2 + a 3 ). 

If the light Higgs h 2 , hi, responsible for the up and down quark masses respectively, arise 
from the same bidoublet h = (1, 2, 2), then they have the same U(l)x charge so that only the 33 
renormalisable Yukawa coupling to h 2 , hi is allowed, and only the 33 element of the associated 
mass matrix will be non-zero. The remaining entries are generated when the U(l)x symmetry 
is broken, via Standard Model singlet fields, which can be either chiral or vector ones. 

In our case, all charge and mass matrices have identical structure under the U(l)x symmetry, 
since all known fermions are accommodated in the same multiplets of the gauge group. The 
charge matrix is of the form 

/ — 2«2 — 4«3 — 3«3 — «2 — 2«3 \ 

-3a 3 2(a 2 - a 3 ) a 2 - a 3 (14) 
\ — a 2 — 2a 3 a 2 — a 3 J 

Then, including the common factor S which multiplies all non-renormalisable entries, the fol- 
lowing hierarchy of dominant effective Yukawa couplings is obtained 

/ £ e |2+6a| £ £ |3a| £ e |l+3a| \ 

A « 5e^ 5e 2 5e , (15) 
V 5el 1+3a l 5 e 1 J 

where we have used a vector-like pair of S'[/(4) <g> SU(2) l ® SU (2)r singlets to Higgs the family 
symmetry. Here e = (< 9 > /M)'" 2- " 3 ' where M is the unification mass scale which governs 
the higher dimension operators, while a unique charge combination a = 0:3/(0:2 — 0^3) appears 
in the exponents of all matrices, as a result of quark-lepton unification. In Eq.|l5| the flavour 
independent suppression factor 5 = (HH) /M 2 results from the operators in Eq.|10[ The factor of 
S actually helps the numerical fit and setting a = 1, e = 0.22 and 5 ~ 0.2 successfully reproduces 
the desired values of the elements of the Yukawa matrix in Eq.[H| Thus, starting with all of 
the "fundamental" Yukawa couplings[| ~ 1, we can reproduce the correct phenomenology of 
charged fermion masses and mixings. 

Let us now look at the case of non-symmetric textures^, with an additional zero in the 23 
entry. It turns out || |n| that it is optimal to consider a scheme in which the dominant operators 



in the Yukawa matrix are O33, 0^ 2 , 0\%, 2 i,0 2 i and Oi 2 , where the last three operators are 
left general and will be specified later. We are also aware from the analysis in ref. that 
O12 must have a zero Clebsch coefficient in the up sector. A combination of two operators 



1 i.e. dimensionless couplings in the full SU(4)®SU(2)L&)SU(2)i{£g>U(l)x gauge theory. 
2 Asymmetric textures have been discussed in a different framework in 021. 
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must then provide a non-zero 2 \ entry to generate a big enough V u b, while an additional much 
more suppressed operator elsewhere in the Yukawa matrix gives the up quark a small mass. At 
Mgut, the Yukawa matrices that we can form from the operators listed above, are of the form 
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(16) 



where only the dominant operators are listed and complex phases have been taken into account. 
The I superscript labels the charge sector and xfj refers to the Clebsch coefficient relevant to 



the charge sector I in the ij position. 0^- are unknown phases and Hy is the magnitude of 
the effective dimensionless Yukawa coupling in the ij position. Any subdominant operators 
that we introduce will be denoted below by a prime and it should be borne in mind that these 
will only affect the up matrix. So far, the known Clebsch coefficients are 
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We have just enough freedom in rotating the phases of -F\,2,3 and -Fi,2,3 to get rid of all but one 
of the phases in Eq.|16|. When the subdominant operator is added, the Yukawa matrices are 
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where we have defined 



Hge^ = H^x^ 21 + H 21 5t ^e** 81 
T = H 21 x^ E e^+H 21 x^ E e 1 ^ 
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(19) 



We may now remove cf) 22 by phase transformations upon i*i,2,3 but (f) 2l may only be removed 
by a phase redefinition of -^1,2,3, which would alter the prediction of the CKM matrix Vckm- 
Thus, 02i is a physical phase, that is it cannot be completely removed by phase rotations upon 
the fields. Once the operators 2 i,0 2 i,Oi 2 have been chosen, the Yukawa matrices at M GUT 
including the phase in the CKM matrix are therefore identified with Hij, H' 22) (p 2l . 
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Table 3: Results of best-fit analysis on models with n = 1 operators. 



This type of non-symmetric texture can be described by a structure of the kind, 

/ 8e Ua 

\= 5e 3 5e lor ' 2 5e M9 I (20) 




where we identify e = A = 0.22 and set 5 ~ 0.1. 

Let us numerically analyse the non-symmetric case we have been discussing, since it goes 
beyond solutions that have been discussed so far in the literature. This is done by performing a 
global fit of possible textures, arising from different operator assignments, to m e ,m M , m u , m c , 
m t , rrid, m s , m,&, ag(Mz), \V u b\, \V c b\ and \V US \ using m T as a constraint. The values of the 8 
parameters introduced in Eq.[L8| {(j) 2l = 0, = H 2 i, H 21 = H 2 i, H 2 2,H 2 2, H l2 , H m ,H 33 ), 
plus as at the GUT scale are determined by the fit. 

The matrices A 7 are diagonalised numerically and \ V u b(M GUT )\, \V US (M GUT )\ are determined 
by Vckm = VulVdl-i where Vjjl, Vdl are the matrices that act upon the (it, c, t) l and (d, s, b) l 
column vectors respectively to transform from the weak eigenstates to the mass eigenstates 
of the quarks. We use the boundary conditions ai(M GUT ) = a 2 (M GUT ) = 0.708, while 
^u,c,t,d, s ,b,e,fi,r, \V US \ and \V u b\ are runQ from M GUT to 170 GeV~ m t using the RGEs for the 
MSSM. The Aj are evolved to their empirically derived running masses using 3 loop QCD®1 



loop QED ||. m e T and A^(m T ) thenQfix tan (3 through the relation |I3| cos/5 = ^^p^ ; where 
v = 246.22 GeV is the VEV of the Standard Model Higgs. There are twelve data points and 
nine parameters so we have three degrees of freedom (dof). The parameters are all varied until 
the global x 2 /dof is minimised. The data used (with la errors quoted) is taken from |Tj|. For 



12 = R , 2 i + 21 = M + A (21) 

we obtain the best fit. The result appears in Table [5| 

3 All rcnormalisation running in this paper is one loop and in the MS scheme. The relevant renormalisation 
group equations (RGEs) are listed in ref. gj. 

4 The superscript e upon masses, mixing angles or diagonal Yukawa couplings denotes an empirically derived 
value, whereas the superscript p denotes the prediction of the model for the particular fit parameters being 
tested. 
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Out of 16 possible models that fit the texture required by Eqs.|T7],|T(], 11 models fit the data 
with x 2 /dof< 3, 5 with x 2 /dof< 2 and 3 with x 2 /dof< 1. The operators listed as Oi2,02i,0 2 i 
describe the structure of the models and the entries H22, #12, H21, cos0, if 33, H22 , H 2 i are 
the GUT scale input parameters of the best fit values of the model. The estimated la deviation 
in as{Mz) from the fits is ±0.003 and the other parameters are constrained to better than 1% 
apart from cos<ft, whose la fit errors often cover the whole possible range. We conclude that 
the x 2 t es t has some discriminatory power, since if all of the models were equally good, we 
would statistically expect to have 11 models with x 2 /dof< 1, 3 models with x 2 /dof= 1 — 2 and 
2 models with x 2 = 2 — 3 out of the 16 tested. 

Summarising, we have combined the idea of a gauged U(l)x family symmetry with that of 
quark-lepton unification within the framework of a string-inspired Pati-Salam model. The non- 
renormalisable operators are composed of a factor {HH) and a factor involving the singlet fields 
9,9 as in Eq.|^. The singlet fields 9,9 break the U(l)x symmetry and provide the horizontal 
family hierarchies while the H,H fields break the SU(4)(8>SU(2)^®SU(2)r symmetry and give 
the vertical splittings arising from group theoretic Clebsch relations between different charge 
sectors. We have studied both symmetric and asymmetric mass textures. In the first case, as 
an example, we have shown how one of the RRR textures can be reproduced by the structure in 
Eq.[L5], where 5 results from the factor (HH)/M 2 and e results from the factor (< 9 > /M)' Q2 ~" 3 ' 
in the operators of Eq.|H]. This results to an accurate reproduction of RRR solution 5 which 
cannot be reproduced by family symmetry alone. The Giudice ansatz WM (RRR solution 3), 
which also cannot normally be achieved, can similarly be reproduced by our scheme. Our 
approach has also been extended to non-symmetric textures, which are motivated by particular 
superstring constructions. For these, we performed a numerical analysis since they are beyond 
the symmetric RRR regime and found that from the possible textures, one of them allows 
a low energy fit to the data, with a total \ 2 of 0.39 for three degrees of freedom. All the 
textures that we present, rely on the existence of Clebsch zeroes, which are a feature of the 
operators in our model. Moreover, the combination of family symmetry and unification has the 
additional advantage that quark and lepton masses are related by Clebsch factors, as in the 
Georgi-Jarlskog scheme, which improves the overall predictive power of the approach. More 
detailed issues, including the superstring construction itself, are deferred to a longer publication 
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